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(1.1) 



Abstract. We consider Hamilton-Jacobi equations, where the Hamiltonian depends discon- 
tinuously on both the spatial and temporal location. Our main results are the existence and 
well— posedness of a viscosity solution to the Cauchy problem. We define a viscosity solution by 
treating the discontinuities in the coefficients analogously to "internal boundaries" . By defining 
an appropriate penalization function, we prove that viscosity solutions are unique. The existence 
of viscosity solutions is established by showing that a sequence of front tracking approximations 
is compact in L°°, and that the limits are viscosity solutions. 



1. Introduction 
In this paper we study the initial value problem 

ju t + H(u x ,a(x),g(t)) =0 t > 0, 



u(x, 0) = uq(x) t = 0. 



where H : M 3 -> R. The functions a : R -» M and g : R+ -> R are called "coefficients", and 
are allowed to be discontinuous. The classical theory for viscosity solutions of Hamilton-Jacobi 
equations, see e.g. jU does not cover the case where the Hamiltonian H is discontinuous. This is 
because the straightforward method of comparing sub- and supcrsolutions does not work if H is 
discontinuous in x or t. 

One application of Hamilton-Jacobi equations with discontinuous coefficients is the optical 
shape-from-shading problem 

(1.2) , 1 =I(x,y), 

1 1 + ul + ul 



where I denotes the intensity of the reflected light, and u = u(x,y) the height of the underlying 
surface. If the gradient of u is discontinuous, then the intensity of the reflected light will vary 
discontinuously. Another related application is the synthetic radar shape-from-shading equation 

(1.3) i =I(x,y). 

\l l + u l + u l 

Also in this case, if the gradient of u is discontinuous, so is /. Note that l|1.3fl can be reformulated 
as the evolution equation 




P(x,y) 



This is reminiscent of (jl . L|) . but unless I(x,y) = f(a(y),g(x)) for some function /, this equation 
is not actually of the type considered in this paper. 

A number of authors have considered Hamilton-Jacobi equations with coefficients that have 
some kind of singularity. In 0] Ishii and Ramaswami considered a boundary, and an initial 
value problem. In this paper |5] the notion of viscosity solution was extended to some classes of 
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discontinuous Hamiltonians. Using the definitions from 0, in Tourin |21j showed that for this 
class of discontinuous Hamiltonians, a unique viscosity solution exists. 

Also, if the Hamiltonian H satisfies some structure conditions, essentially amounting to discon- 
tinuous jumps "in one direction" , there exists a unique solution to the initial value problem, see 
Capuzzo Dolcetta and Perthame or Stromberg 

For shape-from-shading problems, Ostrov ^3 E] showed that provided the viscosity solution 
was bounded and absolutely continuous, then a sequence of approximations, defined by smoothing 
the coefficients, converged to the unique viscosity solution. This viscosity solution was defined 
as the solution of an auxiliary control problem. The technique of using the associated control 
problem (for convex Hamiltonians) was also exploited by Dal Maso and Frankowska in 0E] when 
studying some Hamilton-Jacobi equations with a discontinuous Hamiltonian. 

The Hamilton-Jacobi equation (|1 . 1|> is formally equivalent to the conservation law 



where p = u x . This is an example of a conservation law with discontinuous coefficients, and such 
equations have been extensively studied by a number of authors, see e.g. Klingenberg and Risebro 
[131 EH Gimse and Risebro [B], Klaussen and Risebro In H3 it was shown that smoothing 
the coefficients for the conservation law produced a convergent sequence of solutions, such that the 
limit was a weak solution of the conservation law. For a class of flux functions, Towers |22| proved 
uniqueness within the class of piecewise smooth solutions by using a variant of the Kruzkov 
approach, as well as convergence of monotone difference methods |2Sj- Finally, Seguin and Vovelle 
[T8"| studied a special case of the purely hyperbolic version of (|1 with the flux fuction taking 
the form H(a,p) = ag{p) The authors proved uniqueness of L°° entropy solutions by the Kruzkov 
method |15|. Recently, a quite general theory for conservation laws with discontinuous fluxes 
was established by Karlsen, Risebro and Towers in see also Coclite and Risebro [3] in which 
well-posedness is established for some conservation laws of the type (|1 . 

In this paper we establish both uniqueness and existence of viscosity solutions to The 
uniqueness is demonstrated by using a standard "doubling of the variables" approach, and choos- 
ing the penalization function so that the discontinuities does not influence the technique. This is 
motivated by the techniques used in although we could possibly avoided some of the tech- 
nicalities by using a related approach used by Benth, Karlsen and Reikvam in pQ, for a control 
problem. 

The existence is established by showing that a front tracking scheme, designed for (I1.4JI is well- 
defined and produces a sequence that converges to an appropriate entropy solution. The integrals 
of the approximate solutions to the conservation law, are then shown to be approximate solutions 
to the Hamilton-Jacobi equation, and their limit is a viscosity solution. This program was carried 
out for conservation laws/Hamilton- J acobi equations without x or i dependence by Karlsen and 
Risebro in [TH] , 

To be specific, we make the assumptions: 

(A.l) The Hamiltonian H satisfies the following Lipschitz conditions 



(1.4) 



Pt + H(p,a,g) x = 0, 



(1.5) 
(1.6) 



\H(p, a, g) - H{p, a', g')\ < C (\g - g'\ + \a - a'\) (1 + \p\) , 
\H(p,a ) g)-H(q ) a,g)\<C\p-q\, 



for all p, q, a, a', g, g' 



'', and some positive constant C . 



(A.2) 



(1.7) 




(A.3) 



(1.8) 



-r— (0,a,<?) = and 



d 2 H 



(0,a,g) < 0, 



dp 2 



for all a and g. 
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(A. 4) The map p i— > H(p, a, g) is even, and strictly increasing in the interval (— oo, 0), and strictly 

decreasing in the interval (0, oo), for all a and g. 
(A. 5) We have that 

(1.9) lim BiEl^A > C >0, 

p— >oo p 

for some constant C that is independent of a and g. 

(A. 6) The coefficient a{x) is piecewise continuously diffcrcntiable, with finitely many jump dis- 
continuities in a and a', located at the points x\ < ■ ■ ■ < xm ■ 

(A. 7) The coefficient g(t) is piecewise continuous, with finitely many jump discontinuities in g 
and g', located at the points t\ < t 2 < ■ ■ ■ < ij\r- 

(A. 8) The coefficients a and g are of bounded variation, i.e., 

(l-!0) lfflw(R+) < °°> Mbv(r)<°° 

The coefficients have bounded derivatives away from discontinuities, i.e., 

(1.11) max \a (x)\ < oo, max \g (t)\ < oo. 

x^{x 1 ,...,x AI } t^{ti,...,ijv} 

A viscosity solution of the Cauchy problem (|1 . 1|) is defined as follows: 

Definition 1.1. Fix T > 0, and let u : Ibr = M x [0, T) — > M be a bounded and uniformly 
continuous (BUC) function. We call u a viscosity subsolution of <|l.lf> ifu{x, •) = uq, and if for 
all C 1 functions ip(x,t) such that u — ip has a local maximum at (xo,to) G Ht> then 

(1.12) ip t (x , t ) + min {H (ip x (x 0) t ) , a (xq ) , g (t^)) , H (tp x (x , t ) , a (x£) , g (to))} < 0. 

Analogously, we call u a viscosity supersolution if u(x, ■) = uo, and if for every C 1 function ip such 
that u — ip has a local minimum at (xo,to), then 

(1.13) (ft (x , t ) + max {H (ip x (x ,t Q ) , a (xp ) , g (t^ )) , H ((p x (x , to) , a (x£) , g (to))} > 0. 

Finally, we say that u is a viscosity solution of (|1 . 1|) if it is both a viscosity subsolution and a 
supersolution. 



Setting p — u Xl it formally follows that p solves the Cauchy problem 
(1.14) 



p t + H(p,a(x),g(t)) x = t > 0, 

p(x,0)=p o (x) t = 0, 



where H, a and g are as before. To define an entropy weak solution of (f 1 . 1 4|> we shall need the 
Temple singular mapping ^ defined by 

,,,,, ,w \ • i ^ H(p,a,g) - H(0,a,g) 
(l-!5) *(p,a,9) =sign(p) ^ : . 

By an entropy solution to (|1.14|l we shall mean a function p satifying the following definition: 

Definition 1.2. Let p : Tlx —* K be a measurable function. We say that p is an entropy solution 
of <|1.14[l if the following hold: 

(D.l) p e L 1 (n T )nL°°(n T ), V(p(-,t),a,g(t)) e BV(R) for allt G [0,T). Furthermore, the map 

[0,T) Bin p(t, ■) is i 1 (R) Lipschitz continuous. 
(D.2) The function p is a weak solution of l|1.14|) . i.e., 



(1.16) J J pipt + H(p, a, g)ip x dtdx + j ip(x, 0)po dx = 0, 

for all test functions ip G Co(Ht)- 
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(D.3) The following (entropy) inequality holds for all constants c and all non-negative test func- 
tions if, 

x m +i T 



(1.17) 



J J \p— c\ip t +F(p,x,t,c)ip x dtdx - ^2 J J sign (p- c) H a (c,a(x),g(t))a'(x)ipdtdx 

M T , 

J \H{c,a(x+),g(t))-H(c,a(x-),g(t))\dt>0, 



m=l* 

where we have set xq — — oo, xm+i — oo, and F is given by 
F(p,x,t,c)=sign(p-c)[H(p,a(x),g(t))-H(p,a(x),g(t))], t > 0, xe 
The main result of this paper is summarized by 



Main Theorem. Assume that the assumptions (A.l) - (A. 8) hold. Let uq and vq be two func- 
tions in BUC(M), then there exist corresponding viscosity solutions of Ijl.tjl u = u(x,t) and 
v = v(x,t), satisfying the initial conditions 

u(x 7 0) = u Q (x), v(x,0) = v (x), 

furthermore 

\\u(-,t)-v{-,t)\\ LX3{R) < \\u -v \\ Loo(JSL) , 
for all t € [0,T]. Moreover the functions p and q defined by 

du dv 

P = a~' and ? = «~' 
ox ox 

are entropy weak solutions of (|1.14|) . taking the initial values 

p( ,0) = — , and 9(",0) = — . 

Remark 1.3. From and ^I] we know that 

\\p(;t)-q(;t)\\ Lim < bO»0)-gO,0)|| il(K) . 

The remainder of this paper is organized as follows. In Section [21 we show that viscosity 
solutions are unique by establishing a comparison principle. Then in Section [3] we show existence 
of a solution by first showing that a front tracking algorithm for (|1.4f> is well defined. Then we 
show that this automatically yields a front tracking algorithm for the Hamilton-Jacobi equation. 
The convergence of sequences produced by this algorithm(s) finally gives the existence. 

2. Uniqueness and stability 
This section is devoted to the proof of the comparison principle. 

Theorem 2.1 (Comparison principle). Assume that \(A.l)\ \(A.6)\ \(A.7)\ and \(A.8)\ hold. 

Let u and v be two viscosity solutions of taking the initial values 

u(x,0) — Uo(x), v(x,0) = v (x), 
where uq and vq are in BUC(R). If 

uo(x) < vo(x), x £ R, 

then 

u(x,t) < v(x,t), for x £l andt € [0,T]. 

An immediate consequence of the comparison principle is that viscosity solutions are unique. 
Without loss of generality we shall assume that 

(2.1) x M = 0, 

which in particular implies that a is continuous in the interval (0,oo). To prove Theorem 12 . 1 1 we 
shall need two lemmas, the first of which is 
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Lemma 2.2. Assume that (A.l)\ (A.6)\\(A.7) and (A. 8) hold. Let u be a viscosity solution of 
(II. 1|) taking the initial data uq G BUC(R). Let ip be a continuously differ •entiable test function 
if = p(x,t). If u — if attains a local maximum at a point (xo,T), where Xq ^ {xi, . . . ,xm}, then 

(2.2) p t (x , T) + H {p> x (x , T) , a(x ), g (T~)) < 0, 
while if u — ip attains a local minimum at (xq,T) then 

(2.3) p t (x , T) + H (p x (x , T) , a{x ), g (T~)) > 0. 

Proof. Assume that u — pi has a local maximum at (xq , T) . Without loss of generality, we can 
assume that this maximum is strict. Now for e > 0, define 

<Ps(x,t) = <f(x,t) + 



T-t 

For each e, u — ip e attains a local maximum at a point (x £l t e ), such that 

x £ — > xq, and t e — > T, as e — > 0. 



Since xq £ {x±, . . . , xm } by (A. 6) for sufficiently small e, also x £ £ {xi, . . . , Xm}, and t e £ 
{t\, . . . , t^}. Since u is a viscosity solution, 



We also have that 



and thus 



<Ps,t (x e , t e ) + H (ip EtX (x e ,t e ) , a (x £ ) , g (t s )) < 
tp e , x (x,t) = p> x (x,t) and <p e>t (x,t) = <p t (x,t) + 
ft (x £ ,t £ ) + H (p x (x £ ,t £ ) , a (x £ ) , g (t £ )) < 



{T-t, 



(t - ty ' 

< 0. 



Therefore, by (A.l) (A. 6) and (A. 7) l|2.2|l holds. The inequality (|2.3() is proved by the same 
arguments. □ 

The second lemma we shall need to show Theorem 12. II is 



(A.l)\ \(A~F\ \(A.7)\ and \(A.8)\ hold. Let u and v be two viscosity 



Lemma 2.3. Assume that 

solutions of tQ]l in the strip lis, where < S < t\, taking the initial values 

u(x, 0) = u Q (x), v(x,0) = v (x), 
where uq and vq are in BUC(M). If 

(2.4) u (x) < v (x), x e M, 
then 

u(x, t) < v(x, t), for x € M and t € [0, S]. 
Proof. We assume that (|2.4(l does not hold. Then there exists a constant A > such that 

(2.5) sup [u(x,t) - v(x,t) - 2Xt] > 0. 

o<t<s, xeR 

For e > define the function 

$ e (x, t, y, s) = u(x, t) - v(y, s) - X(t + s) - e(x 2 + y 2 ) 



(2.6) 
where 



(t- S y + (\x\-\y\) 2 1(1 



e \ x 



1 



fV (\x\-\y\Y + e 2 



Kx(x,y), 



K = 2 (max |u| + max \v\) , and x{ x iU) = 



f if x < or y < 0, 
otherwise. 



The function $ e has a global maximum at a point (x £ ,t £ ,y £ , s £ ) £ (Ix [0, S]) 2 , and due to the 
presence of the term Kx, 

x £ > and y £ > 0. 
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e 



Next we show that we can find two constants < C\ < C2 such that 

(2.7) % < x e , y £ < -p, Ci£ 3/2 -e< \x e - y £ \ <C 2 e, and \t e - s e \ < C 2 e, 
Since 

2$ e (x £ ,t £ ,y e ,s £ ) > $ £ (x £) t e ,x e ,t e ) + $ e (y £ , s £ ,y £ ,s £ ) , 

we find that 

(t £ - s £ f + (x £ - y £ f e 3 
2 - h 2- -2 < w (a; s ,i e ) - v (y £ ,s £ ) - u (y £ , s £ ) + v (x £ ,t £ ) . 

£ {x e - Vc) + £2 

Since u and v are globally bounded, the terms on the left hand side of the above inequality is also 
bounded, and thus we can find constants < C3 < C4 such that 

C 3 e 3/2 ~e< \x £ -y £ \ < C 4 e, and \t £ ~s £ \<C A e. 

Next we observe that 

(x £ ,t £ ,y £ ,s £ ) > $ £ ( -|^,t e , s E 

= .(^ ! t E )- u (-L, S£ )-A( tE+Ss )-4-^£ 

This implies that 

£ (xl + ye) + l Q» + ^Jf) <u(x £ ,t e )-v{y e> S e )- U (-^= > t^j + V 

-^yfl s — 2 — +4 - e . 

£2 (x £ -y £ ) 2 +e 2 

The right hand side of this inequality is bounded independently of e, and hence we can find 
constants < C5 < C§ such that 

C5 Ce 
-7= < Ze, y e < 

This finishes the proof of (|2.7() . We can extract some more information from this, since 

-(---)=(y 2 ~x 2 )-^ .J- - e . 

e \ir 2 yf / v e e/ ex 2 ey 2 

Therefore 

(2.8) lim - ( \ - \\ = 0. 
Next, we aim to show the three following limits 

lim — = 0, 

,. x £ -y £ 

(2.9) lim— ^— = 0, 

lim k = 0. 

—0 (x e - y £ f + e 2 

To show this we utilize the inequality 

$ £ {x £ ,t £ ,y £ , s £ ) > $ e (x £ ,t £ ,x £ ,t £ ) 



2 1 

u (x £ ,t e ) - v (x e ,t e ) - 2Xt £ - 2ex 2 - e, 
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which yields 



(t e - s £ ) 2 + (x e - y e ) 



< v (x e , t £ ) - v (y s , s e ) + A (t £ - s £ ) 



Since the viscosity solution v is uniformly continuous, this is easily seen to imply (|2.9(l . Now we 
define the sets 



, - i X e 3x e \ / t e 3t £ 



u F = 



te Jte 

4' 4 



I Ve 3yA I s £ 3s £ 



V B = 



'ys 7y e " 




S e 7s £ ~ 


.4' 4 _ 


X 


4' 4 



Let Xi,e and X2,e be two smooth functions such that 

\ if{x,t)€U s , 



X2,e(2/,s) 



1 if(y,s)eV £ , 



^0 ii(x,t)£U e , " v *' [0 if {y,s)tV e , 
and < X(i,2),e < 1- Using xi )E and xi,e we define the test functions 

<Pe{x,t) = Xl,e(x,t) (y e ,S B ) + A (t + S e ) + £ (x 2 + y £ ) 

, (i - s e ) 2 + - y e ) 2 1/1 , 1 



(2.10) 



(2.11) 



e \x 2 y £ 



tpe(y, s) = X2,e(y, s) (X e ,t e ) ~ A (t E + s) - £ (x 2 £ + y 2 ) 



{t £ - s) 2 - (x - y £ ) 2 1/11 



£ \x 2 y 2 



{X-Ve) +£ 2 



(x £ -y) +£ 2 



Now, since u — if has a local maximum in (x £ , t £ ) and v — ip has a local minimum in (y £ , s e ), and 
a is continuous in x £ and y £ we have that 



and 



<Pe,t [x e ,t e ) + H ((p e<x (x e ,t e ) , a (x £ ) , g (t £ )) < 0, 



^6,3 (ye, s £ ) + H (<ip e y (y £ , s £ ) , a (y £ ) , g (s £ )) > 0. 



Subtracting these, and using first the Lipschitz continuity of H (c.f. (A.l) I, and then the Lipschitz 
continuity of a and g (c.f. (A.6)| (A. 7) I and the fact that t £ and s £ both are less than S < t\, we 
find 

Ve,t (»e, *e) - ^e,« (j/e, s e ) < H {ip e ,x (x £ ,t £ ) ,a(x £ ,g (t e ))) - H (ip e>y (y £ , s £ ) , a (y £ ) , g (s e )) 

< C(\ip e (x £ ,t £ ) - lp £ , y {Ve,S £ )\ 

+ (\g(te) -g(s s )\ + \a(x e ) - a(y e )\) (\<Pe,x (x £ ,t £ )\ + \lfj £ , y {y £ ,s £ )\) 

< C(\(fe (x £ ,t £ ) - lj)e,y {Ve,S E )\ 

+ Qte - S e \ + \x £ - y £ \) (\ip e ,x (x £ ,t £ )\ + \lp £ , y (y £ ,S £ )\) 
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We also have that 

ife.t (x £ , t E ) = A + 2 te Se , 

A,s (Ve, 8 S ) = -A + 2 ^ ^ 



(p SiX [x £ , t £ ) = 2ex e + 2 ^ T - 2" 



f(x £ - y £ ) 2 +e 2 ) 

Ve,y (ye, «s) = -^Ve + 2 ^ + — - 2- 



(Oe - ye) 

and so 

(2.12) ^ (s e , t e ) - (ye, Se) = 2A, 

2/1 1 

(2-13) |<^ £ (x £ ,i £ ) - t/'e.i/ (ye,s £ )| = 2e(a; e + y £ ) + - — + — 



and 



|</>s,x (#e,*e)| + IV'e.y (Ve, S e )\ < 2e (x £ + Jfe) H -5 



' (%e ~ Ve) < \Xe ~ Ve\ ^_ _^ Q 



ln „ 2/1 1 A e 3 (x £ -y £ ) 

(2.14) +_ _ + _)+ 4 — »li — 

£ Ue Vj [(x £ -y £ ) 2 + e^ 

Now we have that 

1 e 3 / 2 1 £ 3 / 2 

if-cp and ^f- Cp 

so the second term in (|2.13|) tends to zero as e — ► 0, and since x £ y £ < C/y/e, the first term in 
(|2.13(l also tends to zero. Thus we have 

(2.15) lim \ip £ (x £ ,t £ ) - ip £ ^ y (y £ ,s £ )\ = 0. 
Furthermore, since 

({x e -y £ f +e* 

as e — > 0, also the last term <|2.14|) will tend to zero as e — > 0. Recalling that \t £ — s £ \ je and 
\x £ — y £ \ je both tend to zero as e — > 0, it follows that 

(2.16) lim (\t £ - s £ \ + \x £ - y £ \) (\<p e>x (x £ ,t £ )\ + {y £ , s £ )\) = 0, 

Now we have established 

2A < lim CI \tp e (x £ ,t £ ) - ip e v (y £ ,s £ )| 
e-*o \ 

= o, 

which is a contradiction, since A > 0. This concludes the proof of Lemma l2~3l □ 

Proof ( of Theorem \2.1\) . If T < t±, the statement of the theorem is identical to that of Lemma l2~3l 
If T > fi, then by Lemma l2~3l u ( ■. f) < v(-,i) for f G [0, ii]. In the strip we can use u(-,ti) 

and u(-,ii) as initial data, and conclude that u(-, t) < v(-,t) also in the strip [ti, t 2 ]- We can repeat 
this argument in each strip [ti, tj+j.] to prove the theorem for any finite T. □ 

By modifying the proof of Lemma 12.31 it is now straightforward to show stability, namely 
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Theorem 2.4 (Stability). Assume that (A.l) (A. 6) (A. 7) and (A. 8) hold. If u and v are 

two viscosity solutions with initial data uq and vq (in BUC(M)) respectively, then 



(2.17) 



IK-,*) -"(•.*)!! 



< IK - foil 



Proof. The proof of this theorem consists in repeating the arguments in the proof of Lemma 12.31 
but starting by using the penalization 

\2 



$ e (x, t, y, s) = u(x, t) - v(y, s) - X(t + s) - s(x 2 + y 2 ) - 



{t-sf + {\x\-\y\Y 



1/1 



S V X 



{\x\-\y\Y 



Kx(x,y) - IKo - Ml 



instead of l|2.6|) . 



□ 



3. The front tracking scheme 

In this section we show existence of a viscosity solution of <|l.lf) ■ This is done by first considering 
a front tracking scheme, which yields approximate solutions to both l|l.ll) and (|1.14|) . As all front 
tracking schemes, this one is based on the solution of Riemann problems, therefore we start by 
detailing this. 

The Riemann problem for the conservation law is the initial value problem where g is constant 
and a and po take two values, i.e., 

, 31 s \pt + H(p,ai) x = 0, p{x,0)=pi ifx<0, 

|p t + H (p, a r ) x = 0, p(x,Q)=p r if x > 0, 

where pi, r and a^ r are constants. Since g is constant, we have omitted the g dependence of H in 
our notation. The entropy solution to this problem is found by finding two p- values p\ such that 
the (scalar) Riemann problem with a flux function H(p,ai), a left state pi and right state p[ is 
solved by using waves of non-positive speed, and the Riemann problem with flux H(p, a r ) and left 
state p' T and right state p r are solved using waves of non-negative speed only. Since the mapping 
p i ► H (p, a) has a global maximum at p = for all a and is even, these states are found as follows: 

If pi < 0: then p[ is in the set [—pi, oo). 
If pi > 0: then p[ is in the set [0, oo). 
If p r < 0: then p r is in the set (— oo,0]. 
If p r > 0: then p' r is in the set (— oo, — p r ]. 

Furthermore, the Rankine-Hugoniot condition implies that 
(3.2) Htf l ,a l )=Hte.,a r ). 

This is still not enough to give a unique solution, and in [Jj the unique entropy solution is deter- 
mined by the (unique) pair (p'i,p' r ) such that 

\P[-Pr\ 

is minimal. For the flux functions considered in this paper, we can always find a unique solution to 
the Riemann problem in this way. This solution will consist of p-waves, over which a is constant, 
and the discontinuity in a, which we call an a-wave. 

Although the solution of the Riemann problem will in general not be a monotone function of 
x/t, the flux, H(p(x,t),a(x)) will be monotone between the two values H(j>i,ai) and H(p r ,a r ). 
This observation can be used to bound the solution of the Riemann problem. Let G ± (ft, a) be the 
two local inverses of H, i.e., 

G + (H(p,a),a) = \p\ , and G~ (H(p, a), a) = - \p\ . 

Since 

m.m{H (p h ai) ,H (p r ,a r )} < H(p(x,t),a(x)) < max{# (pj, ai) , H (p r , a r )} , 
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we have that 

min {G~ (H (p/,aj) ,a r ) , G~ (H (p r ,a r ) ,a/) \pi\ , - \p r \} 

< p(x,i) < max{G + (H (pi,ai) ,a r ) ,G + (H (p r ,a r ) ,a t ) , \pi\ , \p r \) . 
By the assumption (A. 5) we know that the values 

min{G~ (H (p u a{) ,a r ) , G~ (H (p r ,a r ),ai)} and max{G + (H (jpi,ai) ,a r ) , G + (H (p r ,a r ),ai)} 

are finite and bounded. By the special form of H, in particular |(A.3) (A. 2) and (A.4)| we have 
the coarser (but simpler) bound 

(3.4) \p(x,t)\ < G + (mm{H (p h ai) , H (p r , a r )} , max{a ; , a r }) . 

3.1. Front tracking with constant g. We start by defining the front tracking scheme for the 
case where g is constant, this is a variation of the front tracking schemes defined in |13l 
Therefore consider the initial value problem 



p t + H(p, a) x = for x G M, t > 0, 
p(x, 0) = Po(x) for 



(3.5) 
Let 

z(p, a) = — sign (p) (H(p, a) — H(0, a)) and a(a) = H(0,a). 
Since a t— > H(0, a) is non-decreasing, a <— > a(a) is invertiblc. In the (z, a) plane, a waves are 
straight lines of slope ±1. An a-wave connecting two points (zi, ai) and (22, 0^2) have slope 1 if Z\ 
and Z2 are non-positive, and slope —1 if these values are non-negative. If z\ and Z2 have different 
sign, there is no a-wave connecting these points. Since p-waves connect points with the same a 
values, these are horizontal lines in the (z, a) plane. Now fix a (small) number S > 0, and set 
cti = i5, and Zj — j5, for integers i and j. We define pg and a s as piecewise constant functions, 
with a finite number of jump discontininuities, such that 

(3.6) J 10 " ]!™^ '! -*-0. 

bo-P |Li (R) ^o J 

Label the (finite number of) values of p s and a s pi, ■ ■ ■ ,Pm, and Oi, . • • , Ozv respectively. Let otj 
be the jth member of the ordered set 

where m! and M' are chosen such that 

to' < mma(a 5 (x)) < maxct(o 6 (i)) < M'. 

X X 

For ease of notation, set 

a 3 = a' 1 (otj) . 

Next for each otj, we define Zj^ to be the fcth member of the ordered set 



where iV'(j) is such that 

(z_Ar/(j),aj) = -P, and z^ 1 (z N , {j) ,aj) = P, 
where the value P will be determined below. We also set 

Pj,k = z^ 1 (Zjj, a.j) , and //,./ = H {p J:kl cij) . 
Then, for each j, let the approximate flux function H s (p, a) be the piecewise linear interpolant, 

(3.7) H 5 (p, aj) = H jtk + {p- Pj , k ) " JJ -> " lJ . for p E \pj, k ,Pj,k+i}. 

Pj,k+i — Pj,k 

Now the front tracking solution, which we label p s = p 5 {x, t) is constructed as follows. At t = 
we solve the Riemann problems defined by the discontinuities in a s and pg, using the flux function 
H s . The flux function H s is constructed so that p-rarefaction waves are now a series of contact 
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=: P. 



( 3 - 10 ) ... .,. 



discontinuities. Thus the solution of the Ricmann problem is a piecewise constant function of x/t. 
Furthermore, all the intermediate states will be on the grid (pj t ki a j)- At some time t > 0, two 
discontinuities, hereafter called fronts, will collide, giving a new Riemann problem centered at the 
interaction points. We solve this Riemann problem, this gives new fronts that can be tracked until 
the next interaction point and so on. Based on the estimate Ij3.4(l . we also have that 

(3.8) \p 5 (x,t)\ < G + (infH(p 5 (x),a 5 (x)),$upa 5 (x)) 

In order to show that p s is well-defined (for instance we must show that there is no accumulation 
of collision times) we define the Temple functional of a front by 

| A\I/| if w is a it-front, 

(3.9) T(w) = {2 \Aa\ g if w is an a-front, and * r < 

4 |Aa| g if w is an a-front, and > ^i, 

where is defined in (|1 . , and we have included the g value in our notation since we shall need it 
later. For sequence of fronts, define T additively. Next, for the front tracking approximation p 6 , T 
is defined as the sum over all the fronts in p s , and with a slight abuse of notation we write T(p s ). 
From estimates found in |2()l IKSj it follows that t t-y T(p s (-,t)) is non-increasing, which again 
implies that [^(p 5 , a" 5 )] w is bounded. This again implies that p s is well defined, and furthermore 
that there is only a finite number of interactions of fronts for all t > 0. 

Summing up, we have chosen the grid so that the entropy solution to the initial value problem 

p t + H S (p,a s ) x = 0, t>0, xeR 

p(x,0)=p 5 [x), xel, 

can be constructed by front tracking for any time t. Furthermore p will take values that are grid 
points, i.e., for any point (x, f) such that p s and a s are constant at (x, i), 

z (p 5 (x, t) , a s (x)) — Zj_k, for some j and k. 

In particular, this means that 

H s (p s , a 5 ) = H (p s , a 5 ) , almost everywhere. 

For an elaboration and proof of these statements, see . The construction used here differs from 
the construction in in that we have added grid points corresponding to the discretization of 
the initial function po and the coefficient a, instead of choosing discretization that take values on 
the fixed grid in the (z, a) plane. 

3.2. Front tracking in general. Now we can define the front tracking approximation in the case 
where g is not constant, c.f. I|1.14|l . Let g s be a piecewise constant approximation to g, such that 

\\9 S -9\\ L i (R+) -> 0, as <5->0, 
(3.11) ( ) 

\9 \bv((o,t]) - \9\bv({o,t]) ■ 

Define t n such that g is constant on each interval I n = (t n ,t n+1 ]. Assuming that we can define 
front tracking for t < t 11 , we can then use p 1 (• , f") as initial values for a front tracking approximation 
defined in [t n , t n+1 ) . In order to do this we must use a "new" mapping z, since z = z(p, a, g), and 
redefine the grid on which we operate. However, we keep the grid points corresponding to p s (-, t n ). 
In this way, the grid used in the interval I n+1 will contain more points than the one used in I n , 
but since there are only a finite number of intervals I n such that t n < T, for a fixed S, we use a 
finite number of grid points for t < T. If, for t G T n , H s (- 7 -^g 5 ^)) denotes the approximate flux 
function constructed above using H(-,-,g s L n ) and p s (-,t n ), then we have that the front tracking 
construction p s will be an entropy solution of 

p s t +H s (p s ,a s (x),g s (t)) x = 0, t > 0, x€l 

p s (x,0)=p s o (x), xeR. 
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We call the discontinuities in u s fronts, and we have three types, u-fronts, a-fronts and g-fronts 
(that have infinite speed!). 



3.3. Compactness. We aim to show that the sequence {p <5 } i5> q is compact in L 1 , by estimating 



o 5 } 

the variation of ^(p 6 ,a s , g s ). For each time t, such that g s is constant at t, we can view p 5 as 
consisting of a sequence of fronts, u-fronts and a-fronts, and we know that T(p s ) is non-increasing 
in I n . 

For t 6 I n we also know that 

H (p 5 (x, t), a 5 (x),g n ) > inf H (p s (as, t n ) , a s (x), g n ) . 

Now we define the sequence {p n } n>1 by 

(p 1 = G+ (inf, H (p s (x), a^x)^ 1 )) , a, g 1 ) 

\p k = G+ (H (p k -\a,g k )) ,a,g k ) for k > 1, 

where a = inf, a(x), and a = sup, a(x), and <7 5 (i) = g™ for t £ I n . By the continuity of H, 

H(p n ,a,g n+1 )>H(p n ,a,g n )-C\Ag\, 
where Ag = g n+1 — g n . Also 

G + (/i, a ,.9" +1 )<G+(/i,a, 5 ") + C|A 5 |, 
for some constant C. This means that 

<p n + C\Ag\, and thus p n < p 1 + C \g\ BV . 

Now, for t E I n , 

\p 5 (x, t)\<G+ (inf H (p s (x, t n ), a 5 (x), g n ) , a, g n ) 

<G+ (H(F-\a,g n ),a,g n ) 
= p n - 

Hence the sequence of front tracking approximations {p s } is uniformly bounded independently of 
S, i.e., 

(3.13) \p S (x,t)\ <C(=C(p ,H,a,g)). 

where C{- ■ ■ ) is a constant (depending on its arguments only). 

Now we define a "Glimm type" functional, which we shall show that is nonincrcasing in time, 
and then this will imply that the total variation of ^(p s , a, g) is bounded. Set 

(3-14) Q(t)=T(t)\g 5 (-)\ BV([tT]) , 

where with a slight abuse of notation we wr ite T(t) = T(p 5 (-,t)). With these definitions, we can 
state the following lemma. 

Lemma 3.1. There exists a positive constant C, depending only on H, a and g, such that for all 
t > 0, we have that the "Glimm functional" 

(3.15) G(t)=T(t)+CQ{t), 

is nonincreasing in time. 

Proof. The proof of this lemma is very similar to the proof of the corresponding lemma in 3 , and 
we detail only the differences. 

In each interval 7™, we know from |13) that T is non-increasing, and the lemma holds. To prove 
the lemma we must study interactions between p- fronts and g- fronts, and between p- fronts and 
g-fronts. 

Now the proof of the lemma for the interaction of a p-front and an a-front is identical to the 
proof of the corresponding case in , which means that we only must study the interaction of a 
p-front with a g-front. 
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First note that by (A. 3) and (A.l) and (A. 4) there exists a value P and positive constants cq 
and Cq such that 



(3.16) 



\p\ < P => < c < H pp (p,a,g) < C 
P<\p\<C ^0<co< \H p (p,a,g)\ < C Q 



if a < a < a and g < g < <?, 



where g — inf g and g — supg. Now we consider the interaction of a single it-wave and a single 
g-wave. The situation is depicted in Figure ^ For this interaction we claim that 




Figure 1 . The states used in an interaction between a p-wave and a g-w&ve 



(3.17) 



\V(p r ,a,g + ) -V( Pl ,a,g + )\-\V(p r ,a,g )-^(pi,a,g )| 

< C\g+ -g-\ \*(p r ,a,g-) - * (pi, a, g' 



We start by noting that 4>(0, •, •) = * p (0, •, •) = 0, and 

|* (p r , a, g+) - * (p u a,g + ) \ - |* (p r , a, g~) - * (pi, a, g~) \ 

< |tf (p r ,a,g + ) (pi,a,g + ) - * (p r ,a,g-) - * {p U a,g~ 

To prove (|3.17(l we consider different cases. 
Case 1: \pi\ < P and \p r \ < P. Now 

( Pr , a, g+) -* (p,, a, <?+) - * (pr, «, 0~) - * (pi, a, <?") 

Pr 

= f (%(Z,a,g+) -* P (e,a,ff")) ^ 



pi 

Pr 



/" (* p (£, a, 5 + ) " *p (0, a, 5 + ) - *p (£, a, g ) + * P (0, a, <T)) d£ 
(*pp (v, a -,9 + ) - *pp (v, a ,9~)) d-qd£, 



pi 

Pr S 



Pi 

Pr £ g + 



^ppg(r),a,g)dgdrid£, 



pi g- 
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and 



Pr 



pi 

P-r 



(V p (£,a,g )-V p (0,a,g )) d£ 



Pr e 



Pi 



(3.18) 



Moreover observe that 



^ pp (p,a,g) = -sign(p) 



H pp (p,a,g) 
H(0,a,g) 



and 



(3.19) ^ PPg (p, a, g) = sign (p) 



H P p(p,a,g)Hg(0,a,g) - Hp Pg (p,a,g)H(0,a,g) 



To fix ideas, we assume that pi < p r , so that by (A. 2) ^(pi, •, •) < ^{p r , ■, •)■ 
If < pi < p r then 



|*(pr,0,5 + ) -*(pi,a,ff + ) -*(p r ,a,<? )-*(pi,o,5 )| 

<Ci\g+-g- 

= Ci | S + - <r 



drjd^ 



Pi o 



(3.20) 



and 



Pi o 

for some positive constants c% and Ci. Therefore 

\*( Pr ,a,g + ) -*(pi,a,g + ) - * (p r , a, g~) - * (p,, a, cT) | 

< — I* (Pr,a,g~) - * (pi,a,3~)| 

Ci 1 ' v ' 1 

and thus the claim holds. Next, if pi < p r < 0, then 

|*(p r ,a, 5 +) -y(p h a,g + )-V(p r ,a,g-)-y(pi,a,g-)\ 

Pr 



< C 2 \g+ - <?" 



dr)d£, 



pi i 



and 



Pt o 



*(Pr,a,ff )-y(pi,a,g ) > c 2 / / dr?eZ£, 
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for some positive constants C2 and C2 . Hence H3.2()[) holds also in this case. If pi < < p r 
then we write 

\*( Pr ,a,g + ) -* {pi,a,g + ) -V(p r ,a,g-) -*(pi,a,g-)\ 

/ 

<C 3 \g + ~g-\ 



and 



\pi e 
/ 







*(Pr,a,ff )-*(pz,a,0 ) > c 3 



dijd£_ , 



for some positive constants C3 and C3, and l|3.20|) follows. If p r < pi we can use the same 
arguments. 

Case 2: \pi\ < P and |p r | < P. We start by observing that 

Pr- g + 

V(p r ,a,g + )-y(p h a,g+)-y(p r ,a,g-)-$(pi,a,g-)= J J * P3 (£, a, 6>) d<9d£, 



Pi g- 



and 



Pr 

y(p r ,a,g-)-V(pi,a,g-)= J %(jt,a,g~)d^. 



Pi 



Since ^> pg is bounded, we have that 



I* (p r , a, g+) - * (pi, a, g + ) - * (p r , a,g )-^( Pl ,a,g ) | < d \g + -g \ \p r - pi\ , 

and since ^f p > C4 for p $ (—P,P), 

|* (Pr,a,g~) - * (pi,a,g~)\ > c 4 |p r -pj| , 

for some positive constants C4 and C4. Thus (|3.2Q(I follows. 
Case 3: \pi\ < P < p r . In this case we start by writing 

V(p r ,a,g + ) -y(pi,a,g + )-*( Pr ,a,g-) +*(p h a,g-) 

= * ( Pr , a, g+) - (P, a, g+) + * (P, a, g+) - (p h a, g+) 

-■&(p ri a,g-)+*(P,a,g-) - * (P,a,g-) + * (pi,a,g-) 

P « 9 4 



Pi g- 



and 



*(pr,a,$ )-*(pj,a,5 )=*(p r ,o,5 )-*(P,a,g )+tf(P,a,s )-V( P i,a,g ) 



p € 



[ V P (Z,a,g ) d£ + j J * pp (77,a, 5 ) djjdf. 



Pi 



Since the derivatives of \& are bounded, 



|*(p r ,a, fl + ) -*(pj,a,5 + ) -*(p r ,a, S )+tf(p z ,a, fl ) | < C 5 (p r - P + \ 
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for some positive constant C5. Since ^ pp is strictly positive inside [— P, P] and is larger 
than some fixed constant outside this interval, 

P*-p\ 



|*(pr,0,5 )-*(pi,a,fl )|>C 5 p r -P + 
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for some positive constant C5, and thus (|3.2()|1 holds. If pi < 0, then 

/ p2 , 2 

|*(pr,o,5 + ) -*(pi,a,0 + ) -*(pr,o,5 )+*f>i,o,fl ) | < C 8 p r - P + — ^ 



and 

/ p2 + 2 

|*(p r ,a,.g )-*(pi,o,5 )|>c 6 fp r -P + — L 

so <|3.2(Jll holds again. 
Case 4: < P and p r < —P. This is analogous to case 3. 
Case 5: \p r \ < P < pi. This is analogous to case 3. 
Case 6: \p r \ < P and pi < —P. This is analogous to case 3. 
Case 7: pi < — P and p r > P. Now we write 

* (p r , a, <?+) (p u a, g + ) - * (p r , a, g~) + * (p,, a, g") 

= * (jv, a, - * (P, a, + * (P, a, ,g+) - * (-P, a, ,g+) 

+ *(-P,a,ff+) -$( Pl ,a,g + ) -V( Pr ,a,g-) +^(P,a,g~) 

-y(P,a,g-)+y(-P,a,g-)-y(-P,a,g-)+y( Pl ,a,g-) 

Pr 9+ PC 9+ -P9+ 

= // ®P9& a > e ) dd( % + J J J * PPg {v,a,0)d6dr)dZ+ J J 9pg(£,a,6) d9d£, 

p g - -p g - pi gr- 

and 

0<*(p r ,o,fl _ )-« r (pi,o,fl _ ) 
= (j? r , a, - # (P, a, + * (P, a, - vf (-P, a, 
+ *(-Pa,.g-) -*(pi,a,g-) 

Pr P £ -P 

= J * p (^o,fl _ ) d£ . + J j Vpp(v,a,9~) + J * p (f,a,<r) d£. 
p -p p( 

As in the earlier cases, H3.20II is straightforward to show from this. 
Case 8: p r < —P and pi > P. This is analogous to case 7. 
Now the proof of (|3.17(l and thereby of Lemma f3. II is finished. □ 

Let T n —T I and g n — g 5 | /n . Since T is non-increasing in each interval I n , from Lemma |3~T1 
we have that 

T n+1 <T n (l + C\ 

By the Gronwall inequality it follows that 

T(t) <!*(()+) exp k" 
<limT( S )exp(| 3 | w ) 

(3.21) < (|* bo, a, 5(0))| w + 4 |a| w | ff (0)|) e^Uv. 

where the sum in the first line above is over those n such that t n < t. 

This clearly implies that the total variation ^ (p 5 , a s , g 5 (t)) is bounded independently of S and 
t. In particular, this means that the front tracking construction is well-defined, and we have a 
finite number of fronts and interaction of fronts, for < t < T. For a proof of this, see e.g. |13| . 
Furthermore, since \p s \ < C, c.f. (|3.13(l . 

*{C,a,g) <y(u 5 (x,t),a s (x),g s (t)) <*(C,a,g). 

By Helly's theorem, for each fixed t G [0,T], 

\& (p s (-,t), a 5 , g 5 (t)) — ► ip, almost everywhere as S J. 0, 



f +1 ~ 9* 



9^1 
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and by the Lebesgue's dominated convergence theorem also in L 1 (R). Furthermore, by a diagonal 
argument, we can achieve this convergence for a dense countable set {i 7 } C [0, T]. For t 1 in this 
set, define 

p(-,F) = *- 1 (^, a , 5 (^)). 
Hence also p s (-, £ 7 ) converges to some p(-, t 1 ). For any t £ [0, T] we have that 

H^(-^)-/ 2 (-,t)IL 1(K) <||/ l (-,o-/ i (^)|| L1(R) 

+ ||p* (-, f) - p 5 > n\\ L1(R} + \\p s > (., V) - p s * (;t)\\ L1(m , 

where i 7 is such that //(-ji 7 ) — > p(-,t J ). It is easy to show, as in [3], that the map i > u s (-,t) 
is L 1 Lipschitz continuous, so the first and third terms above can be made arbitrarily small by 
choosing Si and 82 small, and the middle term can be made small by choosing t 7 close to t. Hence 
we have that p 5 converges to some function p in L 1 (K x [0,T]). 

Now since p 1— > H(p,a,g) has a unique maximum for p — for all a and we can use the 
same arguments as in PJ to show that p is an entropy solution to l|1.14|l in the sense of (|1.17|) . 
Furthermore, using arguments from and |3J the entropy solution is unique. Summing up, we 
have proved: 



Theorem 3.2. Assume that (ATT) - (A. 8) all hold. Let uq € L 1 (R) be such that the total 
variation of ^(uq, a, g) is bounded. Then there exists a unique entropy solution to (|1.14() p = 
p(x,t). This solution can be constructed as the limit of the front tracking scheme outlined above. 

3.4. Front tracking for the Hamilton Jacobi equation. Now we show how the front track- 
ing approximation to the entropy solution of the conservation law also yields a front tracking 
approximation to the viscosity solution of the Hamilton-Jacobi equation. We start by studying 
the Riemann problem. 

Lemma 3.3. Assume that (A.5)\ (A.4)\ (A. 3) and (A.l) all hold, then the Riemann problem 
for the Hamilton-Jacobi equation 



{u t + H (u x ,ai) = 
u t +H(u x ,a r )=0 
u(a?,0) = uo(0) 
has a unique viscosity solution given by 
(3.23) u{x,t) = u (0) + xp(x 1 t) - t 



if x < and t > , 
if x > and t > , 




x < 
x > 



for t > 0, 



H(p(x,t),ai) 
H(p(x,t),a r ) 

where p(x, t) is the unique entropy solution of the Riemann problem for the conservation law l|3.1|l . 
Remark 3.4. From 13.23JI and the fact that p solves a conservation law, it follows that 

(3.24) u[x,t)=u o (0)-tH o + [ p(x,t)dx, 

Jo 

where Hq — H(p[, a{) — H(p' r , a r ), is an alternative formula for u. This can be shown by observing 
that p = p(x/t), and differentiating (|3.23|) with respect to x, using that pt + H(p, a) x = 0. 

Proof. Let u be defined by H3.23JI . If p is continuous at (x, t) then so is u, and if p has a discontinuity 
moving with speed a, set x = at. By the Rankine-Hugoniot condition, if a ^ 0, 

u (x~ , t) — ito (0) + tap (x~ ,t) — tH (p(x~ , t) , et/ jr ) 

= Uq (0) + t {op (x~ . t) - H {p(x~ , t) , <Z/ iT .) ) 

= u (Q) + t (crp (x + ,t) - H (p(x+ ,t), ai , r )) =u(x+,t) , 

where we use a/ if a < and a r if a > 0. If a = 0, then H(p(0~, t), a/) = H(p(Q + ,t),a r ) which 
gives 

u(0 _ ,t) =u (0)+tH(p(0-,t) ,ai) = u (0) + tH (p (0 + , t) , a r ) =u(0+,r). 
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q t + H (g, a r ) x = 0, q(x,0) 



Thus, u is uniformly continuous. Now let ip be a test function, and assume that u — tp has a 
maximum at (a;o,to)- We proceed by studying two cases. 

Case 1: xq ^ 0. First assume that xq > 0, and let q be a solution of the (scalar) Riemann 
problem, 

fp(to,0+) x<0, 

\p r X > 0. 

Clearly p(x, t) = q(x, t) for x > and t > 0, So by [23 Proposition 2.3] 

V* (x , to) + H (<p x (x , t ) , Or) < 0. 

If Xq < 0, then we replace q by the solution of the Riemann problem with initial data 
given by pi and p(0~ , to) and flux H(q, a{), and reach the same conclusion. 
Case 2: £ = 0. Now by $T23$ it follows that 

p'l = ]imu x (x,t ), and p' r = ]xmu x (x,t ), 

where p\ are the states adjacent to x — of the Riemann solution p. Since u — ip has a 
local maximum, 

Pr < <Px (x ,t Q ) < p[. 
Now since p is a Riemann solution, either 

< p' r < p\ or p' r < p\ < 0. 



We first assume that < p' r < p[, then by (A. 3] 



H (p[, Or) < H (cp x (x ,t ) ,a r ) < H {jp' r , Or) 

(3.25) , / 

= H (Pi,ai) < H (<p x (x , to) ,a t ) < H (p r , aj) . 

Also, since u — ip has a maximum at (0, to): 

<p(0,to) - <p(0,t) < u(0,tp) - u(0,t) 
to — t to — t 

for t < to and to — t sufficiently small. Now by definition 

u(0, to) = u(t, 0) - (t ~t)H (p' h ai) = u(t, 0) - (t ~ t) H (p' r , a r ) , 
and therefore 

rr ( , , „(. S ^ U(0,tp) - U(0,t) ip(0,t O ) - ip(0,t) 

-H (p r , Or) = -H (pi,ai) > j— > j— . 

Passing to the limit as t f to, 

(3.26) -H {p' r , Or) = -if (pj, o,) > (0, to) , 
so 

tpt (0, t ) + H (cp x (0, t ) , Or) < ~H (p' r , a r ) + fT (p' r , Or) = 0. 

By (|3.25|l min{i?(<^ a ;(0, to), o;), H(cp x (0, to), Or)} = H((p x (0, to), a r ) and u is a viscosity 
subsolution. 

Now assume that p' r < p\ < 0, then we find 

H (p' r ,ai) < H (p x (x ,to) ,ai) < H (p[,ai) 

(3.27 ) 

= H (p'r, O r ) < (^i (a!Q, to) , O r ) < H (Pi, Or) . 

Then by (|3~^6l) 

^ (0,to) + H(<p x (0,to),a t ) < -H(p' b ai)+H(p' l ,ai) = 0, 

and in this case, by l|H.27|l. mm{H(ip x (0,t o ), ai), H(ip x (0,t o ), a r )} = H(tp x (0,t ),ai) so u 
is a subsolution. 

To verify that u is also a supersolution, we can use the same arguments. □ 
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In order to define front tracking for Hamilton-Jacobi equations, we now observe that the same 
formula, i|3.23|) holds if H is replaced by the piecewise linear interpolant H 5 . Let a front in p s 
start at a point (xo,to), and let its position be 

x(t) = x + a(t-t ), 

where a is the speed of the front. Now define u 5 (x(t), t) by 

u s {x(t),t) = u 5 (so, to) + (x- x Q ) p\x{t) ± ,t) - (t - t ) H s (/(x(t) ± , a 5 , g\t)) 

{ ' , '" Sl = u s (x ,to) + (t- to) (^(a;(t) ± ,t) - tf 5 (/(a:(t) ± ,a 6 ,/(t))) . 

Between fronts we define u s by linear interpolation. By Lemma l3.3l it follows that u s is the unique 
viscosity solution to the initial value problem 



(3.29) u s t +H 5 (u 5 x ,a s {x),g 5 (t)) = 0, u s (x, 0) = «g(0) + / p'foOjdx, 

Jo 

where the piecewise linear interpolant H s is defined via the front tracking procedure, c.f., Ij3.7|) . 
and Uq(0) is a constant. Furthermore 

(3.30) u s (x, t) = u s (0) - tH b a + / p 5 (x, t) dx, 

Jo 

where H$ = ff' 5 (p' 5 (0 ± , t), a 6 ^), g 5 (t)). Since p 5 converges to p in L 1 , and the traces a" 5 , g s ) 

exists for almost all t, it follows that u 5 also converges in L°°, and its limit is 

(3.31) u(x,i)=u o (0)-tH + [ p(x,t)dx. 

Jo 

Now assume that u — <p has a local maximum at (xq, to), then since converges in i°° to u u s — ip 
has a local maximum at 

(x S , t S ), 

and — > x and t 5 — > t as <5 — * 0. Since u* 5 is a viscosity (sub)solution 
<Pt (x 5 ,t 5 ) 

+ mm{H s (,••• .r- )V (x s <-) ,g s (t s --)),H s (cp x ,a s (x s -+) ,g s (t*"))} < 0, 

Sending 5 I 0, we find that the limit u is a viscosity sub-solution. To prove that u is also a 
supersolution, we can use analogous arguments. Thus u is a viscosity solution of with the 

initial data 

r 

u(x, 0) = uq(0) + / p(x,0)dx. 



Furthermore, for any fixed y we have that 

u\x,t) =u 6 (y,t) + J p 5 (z,t)dz, 

i i | as <$->0, 

r 

u(x,t) = u(y,t) + / p(z,t)dz, 



which means that u x — p almost everywhere. Thus we have proved existence of a viscosity solution, 
and that this solution can be constructed by front tracking, i.e., 



Theorem 3.5. Assume that (A.l) - (A. 8) all hold. For each u € BUC(M) and all T > 0. there 
exists a viscosity solution u £ BUC(M x [0, T}) of <|1.1[1 . Moreover 

p{x,t) = g^OM)) 



is a weak entropy solution of (|1.14() . with the initial condition 

The viscosity solution u can be constructed by the front tracking scheme outlined above. 
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